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Abstract 

As it is well known, every bipartite 2(g) 2 density matrix can be obtained from Bell decompos- 
able states via local quantum operations and classical communications (LQCC). Using this fact, 
the Lewenstein-Sanpera decomposition of an arbitrary bipartite 2 2 density matrix has been 
obtained through LQCC action upon Lewenstein-Sanpera decomposition of Bell decomposable 
states of 2 2 quantum systems, where the product states introduced by Wootters in [W. K. 
Wootters, Phys. Rev. Lett. 80 2245 (1998)] form the best separable approximation ensemble 
for Bell decomposable states. It is shown that in these systems the average concurrence of the 
Lewenstein-Sanpera decomposition is equal to the concurrence of these states. 

Keywords: Qucintum entanglement, Lewenstein-Sanpera decomposition, Con- 
currence, LQCC, Bell decomposable states 
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1 Introduction 

Perhaps, quantum entanglement is the most non classical features of quantum mechanics [|l|, ||] 
which has recently been attracted much attention although it was discovered many decades ago by 
Einstein and SchrOdinger ||l], . It plays a central role in quantum information theory and provides 
potential resource for quantum communication and information processing ggl]. Entanglement 
is usually arise from quantum correlations between separated subsystems which can not be created 
by local actions on each subsystems. By definition, a bipartite mixed state p is said to be separable 
if it can be expressed as 

p = ^Wi pf^ (g) pf^ , lUj > 0, ^Wi = l, 

i i 

where p[^^ and p[^^ denote density matrices of subsystems 1 and 2 respectively. Otherwise the state 
is entangled. 

The central tasks of quantum information theory is to characterize and quantify entangled 
states. A first attempt in characterization of entangled states has been made by Peres and Horodecki 
family [^, |^. Peres showed that a necessary condition for separability of a two partite system is 
that its partial transposition be positive. Horodeckis have shown that this condition is sufficient 
for separability of composite systems only for dimensions 2 2 and 2 (8) 3. 

There is also an increasing attention in quantifying entanglement, particularly for mixed states 
of a bipartite system, and a number of measures have been proposed ||5|, ^ ^ . Among them the 
entanglement of formation has more importance, since it intends to quantify the resources needed 
to create a given entangled state. 



An interesting description of entanglement is Lewenstein-Sanpera decomposition [11|. Lewen- 



stein and Sanpera in [ p^ ] showed that any two partite density matrix can be represented optimally 
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as a sum of a separable state and an entangled state. They have also shown that for 2-qubit systems 
the decomposition reduces to a mixture of a mixed separable state and an entangled pure state, 
thus all non-separability content of the state is concentrated in the pure entangled state. This leads 
to an unambiguous measure of entanglement for any 2-qubit state as entanglement of pure state 
multiplied by the weight of pure part in the decomposition. 



In the Ref. |11], the numerical method for finding the BSA has been reported. Also in 2 ® 2 



systems some analytical results for special states were found in |12]. In |13| we have been able to 
obtain an analytical expression for L-S decomposition of Bell decomposable (BD) states. We have 
also obtained the optimal decomposition for a particular class of states obtained from BD states 
via some restricted LQCC actions. 

In this paper using the fact that, every bipartite 2 (8) 2 density matrix can be obtained from 
Bell decomposable states via local quantum operations and classical communications (LQCC)]T^, 

17, 18|, we obtain the optimal Lewenstein-Sanpera decomposition of an arbitrary bipart 2 2 
density matrix through general LQCC action upon the optimal Lewenstein-Sanpera decomposition 
of BD states of 2 2 quantum systems, where the product states introduced by Wootters in [W. K. 
Wootters, Phys. Rev. Lett. 80 2245 (1998)] form the best separable approximation ensemble for 
BD states. We also show that in these systems the average concurrence of theLewenstein-Sanpera 
decomposition is equal to the concurrence of these states. 

The paper is organized as follows. In section 2 we give a brief review of Bell decomposable states 
together with their separability properties. The concurrence of these states is evaluated in section 
3, via the method presented by Wootters in ||l^ . In section 4 we obtain L-S decomposition of 
these states. By using product states defined by Wootters in we prove that the decomposition 
is optimal. In section 4 we obtain the optimal decomposition for an arbitrary 2 03 2 states by using 
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a general LQCC action which is the main result of this paper, 
conclusion in section 5. 
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The paper is ended with a brief 



2 Bell decomposable states 

In this section we review Bell decomposable (BD) states and some of their properties. A BD state 
is defined by 

4 4 

/> = E^''I^*)(V'*|, 0<p^<l, ^K = l. (2-1) 

i=l 1=1 

where \^pi) are Bell states given by 

l^i> = ITT) + IU», (2-2) 
IV'2) = |TT)-|U)), (2-3) 

i^3) = in) + iiT)), (2-4) 
IV'4) = in) - lit)). (2-5) 

These states form a four simplex (tetrahedral) with its vertices defined hy pi = 1, p2 = 1, Ps = 1 
andp4 = 1 0. 

A necessary condition for separability of composite quantum systems is presented by Peres . 
He showed that if a state is separable then the matrix obtained from partial transposition must be 
positive. Horodecki family Q have shown that Peres criterion provides sufficient condition only for 
separability of mixed quantum states of dimensions 2 2 and 2 3. This implies that the state 
given in Eq. ( |2-1| ) is separable if and only if the following inequalities are satisfying 

Pi<^, for i = 1,2,3,4. (2-6) 
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In the next sections we consider entangled states for which pi > ^. 
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3 Concurrence 

In this section we first give a brief review of concurrence of mixed states. From the various proposed 
measures of quantification of entanglement, the entanglement of formation has a special position 
which in fact intends to quantify the resources needed to create a given entangled state P] . Wootters 



in [10 1 has shown that for a 2-qubit system entanglement of formation of a mixed state p can be 



defined as 



1 1 



Eip) = Hi^- + -V'L^y (3-7) 
where H{x) = — x In x — (1 — x) In (1 — x) is binary entropy and C (p) , called concurrence, is defined 

by 

C{p) = max{0, Ai - As - A3 - A4}, (3-8) 
where Aj are the non-negative eigenvalues, with Ai being the largest one, of the Hermitian matrix 



^ = V VppVp and 

p = {(Jy <^ cry)p*{ay ay), (3-9) 

where p* is the complex conj ugate of p when it is written in a standard basis such as{|tt) ,|Ti)},{|iT) Jii)} 
and ay represent Pauli matrix in local basis {|T) > li)} • 

In order to obtain the concurrence of BD states we follow the method presented by Wootters in 
[p!0|]. Starting from spectral decomposition for BD states, given in ( |2-lj ), we define subnormalized 
orthogonal eigenvectors \vi) as 

\vi) = VPili^i) , {vi I Vj) =Pi6ij. (3-10) 
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Now, we can define states \xi) by 



\^^)=Y.^iJ\''i)^ i = l,2,3,4, 



(3-11) 



such that 



{xi I Xj) = {UTU'^)ij = \i5, 



(3-12) 



where Tij = {vi \ vj) is a symmetric but not necessarily Hermitian matrix. To construct we use 
the fact that for any symmetric matrix r one can always find a unitary matrix U in such a way 
that Xi are real and non-negative, that is, they are the square roots of eigenvalues of rr* which are 
same as eigenvalues of R. Moreover one can always find U such that Aj appear in decreasing order. 
By using the above protocol we get for the state of p given in Eq. (2-1) 



-pi 





p3 
-p4 



V 



(3-13) 



Now it is easy to evaluate Aj which yields 



Xl=Pl, A2=P2, A3=P3, Xa=P4. 



(3-14) 



Then one can evaluate the concurrence of BD states as 



C = - p2 - - P4 = 2pi - 1. 



(3-15) 
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Finally we introduce the unitary matrix U which is going to be used later 

/ \ 

i 



U 



10 
10 
i 



(3-16) 



4 Lewenstein-Sanpera decomposition 



According to Lewenstein-Sanpera decomposition [^], any 2-qubit density matrix p can be written 



as 



P = ^Psep + (1 - A) 



AG [0,1], 



(4-17) 



where psep is a separable density matrix and |^/^) is a pure entangled state. The Lewenstein-Sanpera 
decomposition of a given density matrix p is not unique and, in general, there is a continuum set 
of L-S decomposition to choose from. The optimal decomposition is, however, unique for which A 
is maximal and 



p = A(°P*Viep ^ + (1 - X^°P'^)\iJ^°P'^){ip^°P'^\ , a(°p*) G [0, 1]. 



(4-18) 



Lewenstein and Sanpera in [11| have shown that any other decomposition of the form p = Xpsep + 
(1 — A) !'(/')( '(/'I with p 7^ necessarily implies that A < A*^"^*-* [11]. One should notice that Eq. 

( [4-18 ) is the required optimal L-S decomposition, that is, A is maximal and ps is the best separable 
approximation (BSA). 

Here in this section we obtain L-S decomposition for BD states. Let us consider entangled state 
p which belongs to entangled region defined by pi > ^. We start by writing p as a, convex sum of 
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pure state \ipi) and separable state ps as 

p = Ap, + (l-A)|V^i)(V'i|. (4-19) 

Expanding separable state ps as ps = Yli=iPi iV'i) (V'il ^^id using Eq. ( |2-1| ) for p we arrive at the 
following results 

P'i = -, Pi = ^^ — r for i = 2,3,4, (4-20) 

2' 2(1 -pi) ' ' ' ^ ^ 

and 

A = 2(l-pi). (4-21) 



In the rest of this section we will prove that the decomposition (4-19) is the optimal one. To do so 
we have to find a decomposition for ps in terms of product states \ea, fa), i-e. 

Ps = ^^a\ea, fa) {ea, fa\ (4-22) 



such that the following conditions are satisfied ||ll| 



i) All Aq, are maximal with respect to pa = A^ \ea, fa) {^a, fa\ + i^ — ^) iV'i) (V'll and projector 

Pa — I Cq. , fa) {^a ; fa \ • 

ii) All pairs (Aq, A/3) are maximal with respect to pa/s = Aa \ea, fa) {ea, fal+^p |e/3, fp) {ep, //3I + 
(1 — A) \ipi) {ipi\ and the pairs of projector {Pa, Pp). 

Then according to |11] ps is BSA and the decomposition given in Eq. ( [4-19 ) is optimal. 



Lewenstein and Sanpera in |11] have shown that A is maximal with respect to p and P = \ip) {ip\ 
iff a) if 1-0) ^ nip) then A = 0, and b) if G n{p) then A = (-01 p-^ |0)"^ > 0. They have also 
shown that apair (Ai, A2) is maximal with respect to /) and a pair of projectors (Pi,P2) ifi^: a) if iV'i)) 
IV'2) do not belong to TZ{p) then Ai = A2 = 0; b) if iV'i) does not belong, while 102) G T^ip) then 
Ai = 0, A2 = (V2I p-^ |^2>"^ c) if |0;i), 1^2) G n{p) and (V'll p'^ IV2) = then A^ = p'^ |Vi>"\ 
i = 1, 2; d) finally, if IV'i) , IV'2) G n{p) and (V'll IV'2) / then 
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(4-23) 



Ai= i{i^2\p-'m-\{^Pi\p-'m\)/D, 

A2= HH p-' - \ {Hp~'\i^2) \)/D, 

where D = {ipil iV'i) (^2! P~^ IV'2) - I (V'll P-i IV'2) ?■ 

Now let us return to show that the decomposition given in Eq. ( [4-19 ) is optimal. Wootters 



in [10 1 has shown that any 2 ® 2 separable density matrix can be expanded in terms of following 
product states 



1^1 



) = i (e''^ + e''^ \X2) + e''^ \x^) + e''^ 



J04 



1^3 



Z4) 



) = i(e^^i - e''^ \X2) + 



ifs 12:3) _e^^4 i^.^; 



'^e'"^ \xi) - e'^'' \X2) - e"' 1x3) + e'"^ 1x4) 

provided that Ai — A2 — A3 — A4 < 0. Now, the zero concurrence is guaranteed by choosing phases 
Qi, z = 1, 2, 3, 4 to satisfy the relation Yl,j=\ e^*^^ Aj = 0. 

Now using the fact that for marginal states p^ (located at the boundary of separable region) 
the eigenvalues Aj satisfy constraint Ai — A2 — A3 — A4 = 0, we can choose the phase factors di as 
^2 = 6*3 = ^4 = ^1 + |. Choosing = we arrive at the following product ensemble for ps 

W) = |(-«yp[|V'i) - «y^lV'2) - iyp^lV's) - IV'4)), 

12:2) = |(-i\M 1^1) -^\MIV'2> + «aM|V'3) + \MIV'4)), 

^ ^ ^ ^ (4-28) 

12:3) = |(-iyp[|V'i) +«^/p^lV'2) -«yp^lV'3) + ypllV'4)), 

\za) = li-i^JVAi^i) + i\Jv'2\'4^2) + i\JVz\i^-i) - \/KIV'4)), 
where p[ are defined in Eq. (|4-2q ). 

Let us consider the set of four product vectors {Izq)} and one entangled state iV'i)- In Ref. [|l^] 

it is shown that the ensemble {|^«)} are linearly independent. Evaluating Wronskian determinant 



J04 



(4-24) 
(4-25) 
(4-26) 
(4-27) 
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of vectors \'ipi) and \za) we get Wa = |. This implies that vector is hnearly independent with 
respect to ah vectors \za)- Also evaluating the Wronskian of three vectors \za) and {zp) we 

get 

Wi2 = = ^P2(l - V2), = W2A = Ip^{1 - 2p'^), Wu = W23 = y^il - 2p'^). (4-29) 



Equations ( |4-29| ) shows that in the cases that p has full rank three vectors \za) , \zp) and 1-01) are 
linearly independent. Now we consider the matrices pa = Aq, |zq,) {za\ + {1 — X) \ ipi) {ipi\- Due to 
independence of \za) and l-^i) we can deduce that the range of pa is two dimensional. Thus after 
restriction to its range and defining their dual basis \za) and "01^, we can expand restricted inverse 
p~^ as p~^ = A'^^\za){za\ + (1 — A)~^|0i )(0i| (see appendix). Using Eq. ( |6-49| ) it is easy to see 
that {za\Pa^\za) = ■ This shows that Aq are maximal with respect to pa and the projector Pa- 

Similarly by considering the matrices pap = Aq, \za) {za\ + A^ {zf^) {z^\ + (1 — A) \tpi) {^pl\ and 
taking into account the independence of three vectors \za), {zp) and {'ipi) we see that rang of paf3 is 
three dimensional, where after restriction to its range and defining dual basis \za), |%) and ipi 
we can write restricted inverse p~^ as = Aq^|zq, ){za\ +A'^^\zi3 )( i/jj + (1 — X)^^\ipi )( 0i|. Then 
it is straightforward to get (e„| p'^^ \za) = A'^, p'^^ |%) = A^^ and (zq| p'^^ |%) = 0. 

This implies that the pairs (Aq,, A^) are maximal with respect to /3q^ and the pairs of projectors 
(PajPfs), hence we can conclude that the decomposition given in Eq. ( 4-19| ) is optimal. 



We now consider cases that p has not full rank. Let Pa = for a 7^ 1. In this case Eq. ( 4-29| ) 
shows that the pairs {l^i) , \za)} and also {\z^) , \z^)} for /?, 7 7^ 1, a are no longer independent with 
respect to |0i). In the former case we express |0i) in terms of j^i), |zq) then matrix pia can be 
written in terms of two basis \zi), \za) and after some calculations we get (^zi\p^^\zi'^ = Asi^ilz^^ 
Za\Pia\za) = ^^^^ and (zilpfil^a) = where Tiq = AiAq + 2(1 - A)(Ai + Aq). By 
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using the above results together with Eqs. ( 4-23| ) we obtain the maximality of pair (Ai, Aq,) with 



respect to pia and the pair of projectors {Pi,Pa). 

Similarly for latter case we express in terms of {zp), \z^) then matrix pp^ can be written 
in terms of two basis \zp), \z^) and we get ^ -2/3 1 1-^/3^ = Aod^iL^^ (-^71^/37^7) = ^'^"'r^^ '^^ ^^'^ 
1/^/37 1 -^7) ~ where F/s^ = A/3A7 + 2(1 — A)(A^ + A7). Again using the above results 



together with Eqs. ( [4-23 ) we obtain the maximality of pairs (A^, A7) with respect to pp^ and the 



pairs of projectors [Pp,P^). 

Finally let us consider cases that rank p is 2. Let Pa = Pp = ^ for a,l3 ^ 1. In this cases we 
have \za) = \zp) and \zi) = \z^) for 7 7^ l,a,(3. It is now sufficient to take \zi) and \za) as product 
ensemble. But Eq. ( 4-29| ) shows that these vectors are not independent any more, so that we can 



express {tpi) in terms of l^i) and \za), therefore, matrix pia can be written in terms of two vectors 
\zi) and \za) and we get after some calculations (zi\p^^\zi\ = Asi^iiz^^ (za\pi^\za ~ 



r 



and (zi\p-^^^\zo,) = ^p'*' where Eiq, = AiA^ + 2(1 — A)(Ai + Aq). Using the above results together 



with Eqs. (4-23) we deduce the maximality of pairs (Ai, Aq,) with respect to pia and the pairs of 



projectors (Pi,Pq,). 

5 Behavior of L-S decomposition under LQCC 

In this section we study the behavior of L-S decomposition under local quantum operations and 

classical communications (LQCC). A general LQCC is defined by |]l5|, ^ 

{A0B)p{A<g,B)'i 

p = J—, (5-30) 

Tr{{A^B)p{A0By) 

where operators A and B can be written as 

A(S)B = Ua f^'"'"" Ub Z'^'^'", (5-31) 
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where Ua and Ub are unitary operators acting on subsystems A and B, respectively and the 
filtration / is defined by 

jM,a,m^ ;u(/2 + am.cj), 

(5-32) 

fu,b,n^ jy(/2 + 6n.(j). 

As it is shown in Refs. |15, 16|, the concurrence of the state p transforms under LQCC of the form 
given in Eq. ( |5-30 ) as 

(5-33) 



Tr{{A(g)B)p{A(S)B)^) 
Performing LQCC on L-S decomposition of BD states we get 



{A0 B)p{A^ By 
Tr{{A(g) B)p{A0 B)'f) 



AV, + (1-A')|V'>(^" 



with Pg and {ip') defined as 



respectively, and A' is 



Ps -- 



_ {A(S)B)p,{A^By 

Tr{{A(S)B)ps{A(S)B)^y 
^ {A0B)\i;i) 

v/('0i|(AAt0Sst)|Vi)' 
Tr{{A0 B)ps{A(S) B)^) ^ 



Tr{{A(g> B)p{A(g) B)^ 
Using Eq. (|^), we get for the weight of entangled part in the decomposition ( ^-34]) 



(1 _ V) = (V'ilM^-Bi3t)|V,,) _ 



(5-34) 

(5-35) 
(5-36) 

(5-37) 



(5-38) 



Tr{{A^ B)p{A^ B)'f) 

Now we can easily evaluate the average concurrence of p' in the L-S decomposition given in ( 5-34|) 

^2 z.2^i_q2)(i_^2) 



(l-A')C(IV')) 



:i-A)c(iV'i», 



(5-39) 



Tr{{A(g> B)p{A(g> B)^) 

where, by comparing the above equation with Eq. ( 5-33| ) we see that (1 — X)C{\i(j)) (the average 
concurrence in the L-S decomposition) transforms like the concurrence under LQCC. 
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Now we would like to show that the decomposition given in Eq. ( 5-34| ) is optimal. To do so, 
we perform LQCC action on matrices pa = Aq \za) {za\ + (1 — A) \ipi) {ipil and get 

{A0 B)pa{A0 B)^ 



Tr{{A0B)pa{A0B)'() 



A'^K){z'^\ + {l-X')\^[){^[\ 



(5-40) 



where 



iA<E)B) \Za 



and 



A' 



'{za\{AA^ (^BB^)\za) 
{AA'f ®BB'f)\za) 



Ar 



(5-41) 



(5-42) 



Tr{{A0 B)pa{A0 B)^) 
Using the fact that LQCC transformations are invertible [16, 17, 18], we can evaluate p'^ ^ as 



p'^ ' = TriiA B)pa.{A 5)t) (^t ® Bt)-ip-i(^ ^ B)-\ 



(5-43) 



Using the above equation and Eq. ( 5-41 ) we get 

Pa \^a) 



, I j-i I ^, , _ Tr{{A ® B)p^{A ® S)t) , 



{zc]{AA\ ®BB\)\zo) 



(.^a \ Poi l^ct) — Aq, 



(5-44) 



Equation ( 5-44 ) shows that A^s are maximal with respect to p'^ and the projector P'^. 
Matrices p^p = A^ \za) {za\ + A^ \zp) {zp\ + (1 — A) 1^"!) (V'll transform under LQCC as 



I 

Paf3 



{A® B)pa,^fi{A® B) 



Tr{{A(®B)p^f^{A®B)^) 



(5-45) 



where 



{A(^B)\z^^l3) 



and 



a: 



{ZaAi^^^^ ®BB\)\z^^p) 

liAA'f ®BB'f)\z^^p) 
Tr{{A®B)p^piA0B)'i) 



Aa,l3- 



(5-46) 



(5-47) 
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We now consider cases that p is full rank. In these cases we have already showed that all vectors 
\za) , \zp) and are linearly independent. Using the above results together with invertibility of 
LQCC actions we arrive at the following results 

1""^ I = (.,|(AAtUBt)|J) <^^l ""^ = <^-^«' 

{za\ (AAl SSt) \za) {z^\ {AA\ ®BB\)\zp) 



Equations. ( |5-48| ) show that the pair (AJ^, A^) are maximal with respect to p'^ ^ and the pair of 
projectors {P^,P'^). For other cases that p is not full rank we saw that there is some dependency 
between three vectors \za) , {zp) and iV'i) such that {za\ p~p) {zp) / 0. This implies that in general 



Pap 



z'n) 7^ 0. In this cases in [^] we have shown that under restricted LQCC actions for 



which A = B, the optimality of the decomposition given in (5-34) will be achieved 



6 Conclusion 

We have derived Lewenstein-Sanpera decomposition for BD states and have showed that for these 
states the average concurrence of the decomposition is equal to their concurrence. It is also shown 



that product states introduced by Wootters in |1C] form BSA ensemble for these states. By per 



forming LQCC action on these states we have been able to obtain optimal decomposition for all 
2 2 systems. It is also shown that for these states the average concurrence of the decomposition 
is equal to their concurrence. 
Appendix 
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Let us consider the set of linearly independent vectors then one can define their dual 
vectors { 4>i^} such that the following relation 

h \<t>j) = Sij (6-49) 



hold. It is straightforward to show that the and their dual { (f>ij} posses the following 

completeness relation 

Y,m{'Pi\=I, ^|0,)(0,|=/. (6-50) 

i i 

Consider an invertible operator M which is expanded in terms of states \4>i) as 

M = J2a^j\(|)^){cP,\ (6-51) 

i 

Then the inverse of M denoted by can be expanded in terms of dual bases as 

i 

where bij = {A~^)ij and Aij = aij. 
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